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. Abstract 

pLn ' In this study we investigate the bound states of the Hamiltonian describing a 

quantum particle hving on three dimensional straight strip of width d. We im- 
pose the Neumann boundary condition on a disc window of radius a and Dirichlet 
boundary conditions on the remained part of the boundary of the strip. We prove 
Ph . that such system exhibits discrete eigenvalues below the essential spectrum for any 

a > 0. We give also a numeric estimation of the number of discrete eigenvalue as 
I a function of — . When a tends to the infinity, the asymptotic of the eigenvalue is 

given. 
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cn 1 Itroduction 
(N 

I The study of quantum waves on quantum waveguide has gained much interest and has 



been intensively studied during the last years for their important physical consequences. 
The main reason is that they represent an interesting physical effect with important 
applications in nanophysical devices, but also in fiat electromagnetic waveguide. See the 
monograph [10] and the references therein. 

Exner et al. have done seminal works in this field. They obtained results in different 
contexts, we quote [21 El [HI E] • Also in [121 [131 [HI research has been conducted in this 
area; the first is about the discrete case and the two others for deals with the random 
quantum waveguide. 

It should be noticed that the spectral properties essentially depends on the geometry 
of the waveguide, in particular, the existence of a bound states induced by curvature 
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[31 El [71 [Sj or by coupling of straight waveguides through windows [H [10] were shown. 
The waveguide with Neumann boundary condition were also investigated in several pa- 
pers [T^l [T7] . A possible next generalization are waveguides with combined Dirichlet and 
Neumann boundary conditions on different parts of the boundary. The presence of differ- 
ent boundary conditions also gives rise to nontrivial properties like the existence of bound 
states. The rest of the paper is organized as follows, in Section 2, we define the model and 
recall some known results. In section 3, we present the main result of this note followed 
by a discussion. Section 4 is devoted for numerical experiments. 

2 The model 

The system we are going to study is given in Fig 1. We consider a Schrodinger particle 
whose motion is confined to a pair of parallel plans of width d. For simplicity, we assume 
that they are placed at 2; = and z = d. We shall denote this configuration space by VL 

VL = R^ X [0,rf]. 

Let 7(a) be a disc of radius a, without loss of generality we assume that the center of 
7(a) is the point (0,0,0); 

7(a) = {(x, 0) e M^; + < a"}. (2.1) 

We set r = 9fi\7(a). We consider Dirichlet boundary condition on F and Neumann 
boundary condition in 7(a). 

2.1 The Hamiltonian 

Let us define the self-adjoint operator on L?'{VL) corresponding to the particle Hamiltonian 
H . This is will be done by the mean of quadratic forms. Precisely, let go be the quadratic 
form 

qoU.g) = I Vf-V~gd^x, with domain Q{q,) = {/ e H\n)- f\V = 0}, (2.2) 
Jn 

where H^{Q) = {/ G L^(f2)|V/ G L^(f2)} is the standard Sobolev space and we denote by 
/[F, the trace of the function / on F. It follows that go is a densely defined, symmetric, 
positive and closed quadratic form. We denote the unique self-adjoint operator associated 
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Dirichlet boundary condtion 



Neumann boundary condition 




Figure 1: The waveguide with a disc window and two different boundaries conditions 

to go by H and its domain by D{Q). It is the hamiltonian describing our system. From 
[18] (page 276), we infer that the domain D{VL) of H is 

df 



and 



D{Q) = {/ G H\Q)- -Af e L'{n), /rr = 0, ^ Ma) = o} 

Hf = -Af, WfeD{n). 



2.2 Some known facts 

Let us start this subsection by recalling that in the particular case when a = 0, we get 
H^, the Dirichlet Laplacian, and a = +oo we get H°°, the Dirichlet- Neumann Laplacian. 
Since 

H = (-Ak2) 01® 10 (-A[o,d]), on L2(r2) ® ^^([0, d]), 
( see [18]) we get that the spectrum of is [(^)^, +oo[. Consequently, we have 
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Using the property that the essential spectra is preserved under compact perturbation, 
we deduce that the essential spectrum of H is 



An immediate consequence is the discrete spectrum lies in (f )^] • 



2.3 Preliminary: Cylindrical coordinates 

Let us notice that the system has a cylindrical symmetry, therefore, it is natural to 
consider the cylindrical coordinates system {r^O^z). Indeed, we have that 

L2(1], dxdydz) = L^QO, +oo[x [0, 27r[x [0, d],rdrdedz), 

We note by Qr, the scaler product in L'^{Q, dxdydz) = 
L^QO, +oo[x [0, 2tt[x [0, d],rdrd9dz) given by 



(/,fi')r = / fgrdrdOdz. 

-']0,+oo[x[0,27r[x[0,(i] 

We denote the gradient in cylindrical coordinates by V,.. While the Laplacian operator 
in cylindrical coordinates is given by 

ld_ d_ (f_ 
r dr dr dO'^ dz'^ 

Therefore, the eigenvalue equation is given by 



A^,^,^ = r-(^-) + -T^ + — • (2-3) 



-/\r,eJ{r.e,z) = Ef{r,e,z). (2.4) 

Since the operator is positive, we set E = k"^. The equation (12 ■4p is solved by separating 
variables and considering f{r,6,z) = ^{r) ■ il){9)x{z). Plugging the last expression in 
equation (12 ■4p and first separate x by putting all the z dependence in one term so that ^ 
can only be constant. The constant is taken as —s^ for convenience. Second, we separate 



r 

axial symmetry and the solution has to be 27r periodic and single value in ^, we obtain 



the term ^ which has all the Q dependance. Using the fact that the problem has an 
me 

%- should be a constant —n^ for n G Z. Finally, we get the following equation for (p 



1 

^"{r) + -^\r) + [e -s^- -^Mr) = 0. (2.5) 

We notice that the equation (12. 5p . is the Bessel equation and its solutions could be ex- 
pressed in terms of Bessel functions. More explicit solutions could be given by considering 
boundary conditions. 



4 



3 The result 



The main result of this note is the following Theorem. 

Theorem 3.1 The operator H has at least one isolated eigenvalue in [(^)^, (f )^] /^'^ 
any a > 0. 

Moreover for a big enough, if \{a) is an eigenvalue of H less then —, then we have. 




Proof. Let us start by proving the first claim of the Theorem. To do so, we define 
the quadratic form Qo, 



Qo(/, g) = (V/, Vg)r = [ {drfdrg + \dgfdeg + djd,g)rdrdedz, (3.7) 

^]0,+oo[x[0,27r[x[0,d] 

with domain 

Voin) = {/ e L'^{VL,rdrdedz)]Vrf G L'^{Vt,rdrdedz)-J\T = O} . 
Consider the functional q defined by 



qm = Qom-{^fmhin,rdrdo<i.y (3.8) 

Since the essential spectrum of H starts at (-7)^, if we construct a trial function $ G 2^o(^) 
such that has a negative value then the task is achieved. Using the quadratic form 
domain, $ must be continuous inside Q but not necessarily smooth. Let x be the first 
transverse mode, i.e. 



/ X I A/jiSmi^z) if zE(0,d) „x 



otherwise. 

For ^{r,9,z) = (p{r)x{z), we compute 

g[$] = {Vr^x,^r^x) -i'^yyxWh 



0,+oo[x[0,27r[x[0,d] 
2 

L'^{[0,+oo[,rdr) 



^) \\'rA\\L'^{n,rdrdBdz)^ 



(lx(^)n^'(r)P + W{r)\W{z)\') rdrdOdz - {^hxllhin 



27r||v?'ii2 
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Now let us consider an interval J = [0, b] for a positive b > a and a function (p G 
5 ([0, +oo[) such that (p{r) = 1 for r G J. We also define a family {(p^ : r > 0} by 

'^"^'■^ " { ^(6 + r(lir-ln6)) if r>b.^ ^^'^^^ 

Let us write 

lbrlU2([o,+oo),rdr) = / |¥?Ur)|Vrfr, 

7(0,oo) 



/ T'^\(p'{b + T{\nr — lnb))frdr, 

J{b,+oo) 

/ —\ip'(b + T{lio.r — lnb))\'^rdr, 

J{b,+oo) ^ 

[ -|v9'(6 + r(lnr-ln6))|2rfr, 

J(b.+oo) 



T 



T 



(6,+oo) 

v'is)\'ds = Ty\\l^^^,^+oo)y (3.11) 



'(0,+oo) 

Let j be a localization function from C^{0, a) and for r, e > we define 

$,,,(r, z) = p}^{r)[x{z) + ej{rf] = Vr{r)x{z) + ^,ef{r) = $i,,,,(r, z) + $2,r,.(r). (3.12) 

g[<l>] = g[<I>i,^,^ + $2,rJ 



Qo[*^'l,r,e] - {-^)'^\\^l,T,e\\h(^Q^rdrdedz) + Qo[^2,T,e] " (^)^ II '^2,T,e || i2(f^ ^.^^^0^^) 



Using the properties of noting that the supports of ip and j are disjoints and taking 
into account equation (13.111) . we get 

TT 

= 2nTy'\\L2^o,+oo)-87Tde\\f\\l2^o^+^-^+26\{2\\jj'\\lL2(^o,o^^^^^ 

(3.13) 

Firstly, we notice that only the first term of the last equation depends on r. Secondly, 
the linear term in e is negative and could be chosen sufficiently small so that it dominates 
over the quadratic one. Fixing this e and then choosing r sufficiently small the right hand 
side of (13.131) is negative. This ends the proof of the first claim. 

The proof of the second claim is based on bracketing argument. Let us split L'^{Q, rdrdOdz) 
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as follows, L'^{Q,rdrd9dz) = L'^{il~ ,rdrd9dz) © L'^{il+ ,rdrd9dz), with 

= {{r,e,z) e [0,a] X [0,27r[x[0,rf]}, 

Therefore 

Here we index by D and depending on the boundary conditions considered on the 
surface r = a. The min-max principle leads to 



Vd' 



Hence if H exhibits a discrete spectrum below — , then H do as well. We mention 

d'^ 

that this is not a necessary condition. If we denote by Xj{H~'^), Xj{H~'^) and Xj{H), 
the j'-th eigenvalue of H~'^, H'''^ and H respectively then, again the minimax principle 
yields the following 

XAHa'"") < < ^AHa'"") (3-14) 

and for 2 > j 

A,-i(i/-'^) < X,{H) < X,iH-'^). (3.15) 



has a sequence of eigenvalues [U |T9], given by 



X 



k,n,l 



2d 

Where Xn.i is the l-th positive zero of Bessel function of order n ( see [H [19]) . The 
condition 



yields that /c = 0, so we get 



Xk,n,i < 35", (3.16) 



\2dJ \ a 

a 

We recall that Xn,i are the positive zeros of the Bessel function J„. So, for any A (a) 
eigenvalue of H, there exists, n, I, n', I' G N, such that 



The proof of (13.61) is completed by observing by that Xn,i and Xn',i' are independent from 
a. In Figure 2, the domain of existence of Xi{H), X2{H) and Xs^H) are represented. □ 

4 Numerical computations 

This section is devoted to some numerical computations. In [11] and [16], the number of 
positive zeros of Bessel functions less than A is estimated by — which is based on the 
approximate formula for the roots of Bessel functions for large / is 

a;„,z~(n + 2/-i)|. (4.18) 

taking into account (I3.14p . we get that for, d and a positives such that — < A* = 1, 9276, 
H has a unique discrete eigenvalue. 



Area of existence of the first three eigenvalues of H for d=7i Threshold of appearance of eigenvalues 




0.2 1 1 1 0^ ' ' ' 1 

5 10 15 2 4 6 8 10 

a a 



Figure 2: We represent a i— (^)^ + (^)^ where x(l), x(2), x(3) are the first three zeros 
of the bessel functions increasingly ordered. 
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Figure 3: The number of the eigenvalues of the operator function oi X = a/d. 
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